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We present an approach to spin dynamics by extending the optical Bloch equations for the driven 
two-level system to derive microscopic expressions for the transverse and longitudinal spin relaxation 
times. This is done for the 6-level system of electron and hole subband states in a semiconductor or 
a semiconductor quantum structure to account for the degrees-of-freedom of the carrier spin and the 
polarization of the exciting light and includes the scattering between carriers and lattice vibrations 
on a microscopic level. For the subsystem of the spin-split electron subbands we treat the electron- 
phonon interaction in second order and derive a set of equations of motion for the 2x2 spin-density 
matrix which describes the electron spin dynamics and contains microscopic expressions for the 
longitudinal (Ti) and the transverse (T2) spin relaxation times. Their meaning will be discussed in 
relation to experimental investigations of these quantities. 



PACS numbers: 72.25.Fe,72.25.Rb,78.47.+p 



I. INTRODUCTION 



The Bloch equations, originally formulated as equa- 
tions of motion (EOM) for magnetic moments have 
turned out to apply in general to the dynamics of quan- 
tum mechanical two- level systems . One prominent ex- 
ample are the optical Bloch equations (OBE) in atomic or 
semiconductor physics with the components of the Bloch 
vector composed of the entries of the density matrix for a 
driven two-level system under excitation by a scalar light 
field (see e.g. Ref. 0). Usually carrier scattering is ac- 
counted for by adding phenomenological damping terms 
connected with a longitudinal (Ti) and a transverse (T2) 
relaxation time. In the context of OBE, T\ character- 
izes the decay time of the population inversion or the 
relaxation into an equilibrium distribution, while T2 is 
the timescale on which the coherence between exciting 
light and optical polarization gets lost. A further evolu- 
tion of the OBE are the semiconductor Bloch equations 
(SBE), which were formulated to describe optical phe- 
nomena in semiconductors under intense excitation by in- 
cluding many-particle terms due to Coulomb interaction 
between the carriers Q. These equations yield a micro- 
scopic formulation of Ti and T2 caused by carrier-carrier 
[j| or carrier-phonon scattering 0, ■ In spite of their 
successful application to carrier dynamics, the OBE and 
SBE, in their original form, are not capable to contribute 
to the current topic of spin dynamics in semiconduc- 
tors. Recently, this shortcoming was partially overcome 
by extending the SBE with respect to the spin degree- 
of-freedom of the carriers (including spin-orbit coupling) 
and the polarization degree-of-freedom of the exciting 
light 0, necessary to create a non-equilibrium spin dis- 
tribution by optical orientation . These extended SBE 
(derived by applying the Hartree-Fock truncation) are 
restricted, however, to the coherent regime and hence 
fall short of describing scattering as origin of spin relax- 
ation and spin decoherence, which are key issues of spin- 
tronics and quantum computation |lfj|. On the other 
hand, the structure of these equations resembles those 



used in the phenomenological approach of spin dynamics 
0, , thus indicating the possibility of arriving at a 
microscopic approach to spin relaxation in the language 
of Bloch equations. 

It is the aim of this paper to provide a microscopic 
formulation of spin dynamics in semiconductors and 
semiconductor heterostructures. We do this by start- 
ing from the extended OBE for the 6-level system of 
electron and hole subband states, containing the spin 
of the carriers and the polarization of the exciting light 
(this corresponds to taking into account only the single- 
particle contributions to the SBE of Ref. 0) and include 
the electron-phonon interaction as a possible scattering 
mechanism. For the subsystem of the conduction band 
states (spin-split by spin-orbit coupling) we formulate 
the full dynamics as a set of EOM for the 2x2 spin- 
density matrix and the phonon assisted density matrices. 
By using a correlation-based truncation scheme in sec- 
ond order Born approximation, we derive the scattering 
rates (in the Boltzmann limit and beyond) to arrive at 
equations describing the electron spin dynamics includ- 
ing relaxation (due to electron-phonon interaction) on a 
strictly microscopic level, while existing theories (see e.g. 
Ref. Il2|) are a mixture of microscopies and phenomenol- 
ogy. We want to stress also that, regarding the creation 
of a non-equilibrium spin population, our theoretical con- 
cept differs from some experimental situations: in our 
OBE a non-equilibrium spin polarization is due to opti- 
cal orientation, while in spintronic devices it is usually 
created by spin injection |l3j|. However, this difference 
will not become relevant in the context of this paper con- 
centrating on the spin relaxation due to carrier-phonon 
interaction. 

This paper is organized as follows: In Sec. [D] we in- 
troduce the system Hamiltonian, formulate the full dy- 
namics of the system without truncation and derive the 
EOM for the electron subsystem. In Sees. Hill and IIVI we 
present the correlation-based truncation scheme used to 
achieve a closed set of equations for the entries in the 
2x2 density matrix related to the spin-split conduction 
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band states. It represents an extension of the coherent 
OBE for the spin-density matrix by contributions due to 
electron-phonon scattering. In Sec. we relate the dy- 
namics of the density matrix with those of experimental 
observables and discuss the meaning of the corresponding 
spin relaxation times. Finally, we draw the conclusions 
of our results and give an outlook. 

II. SPIN-DEPENDENT OBE INCLUDING 
CARRIER-PHONON INTERACTION 

The Hamiltonian of the system is formulated in second 
quantization using the notation of Ref.0. We restrict our 



discussion to the case of a quantum well structure (QW), 
but the equations can be formulated in the same way for 
a bulk semiconductor. We consider a six-level system 
consisting of states from the spin-split lowest electron 
subband (with angular momentum or pseudospin indices 
m c — ±i) and the corresponding heavy (m v = ±|) and 
light hole (m v — subband at wave vector k under 
excitation by a light field of arbitrary polarization and 
due to carrier-phonon interaction 

H = Ho + H-light + Ti-phonon ■ (1) 

In the following, we address the individual contributions 
to H. Following Ref. [H we adopt the diagonal form of 



Wo= £ e < (k') C? t „,(k') C< (k')+ £ e m; (k')^(k')^,(k') 



(2) 




FIG. 1: Sketch of subband splitting due to spin-orbit cou- 
pling. 

written with annihilation operators for electrons (holes) 
c. m ' c (k') (u m ^(k')) and corresponding creation opera- 
tors. The time dependence of these operators is under- 
stood. The single-particle energies are denoted as t m > (k') 
(e m ^(k')) for electrons (holes). Although the structure 
of these single-particle contributions is formally equiv- 



alent to a multisubband approach (see e.g. Ref. X\M 
the physical content differs: In the multiband case the 
eigenenergies describe different bands separated by an 
energy gap (e.g. first and second electron subband). In 
contrast, we deal with subband states whose spin de- 
generacyis lifted due to k-dependent spin-orbit coupling 
ll^LiilLLII caused by bulk inversion asymmetry (BIA or 
Dresselhaus term)[l£| and/or structure inversion asym- 
metry (SIA or Rashba term')[T^. |20| (see Fig. QJ, which 
are associated with spin precession. The diagonal form 
of 7io means that the angular momentum or pseudospin 
is defined with respect to the direction of the wave vector 
k . This particular dependence of the pseudospin orien- 
tation on the direction of k is visualized for the Rashba 
spin-orbit interaction in Fig. [3 In dipole approximation, 
the interaction of the light field with the electrons and 
holes reads 



n light = - J2 { E W • d < mi GO <L> e GO vl, (k') + h.c.} , 



(3) 



where E(£) is the electric field vector of the exciting light 
and d™, m , (k') is the dipole matrix element between the 
two subband states with pseudospin index m' c and m' v . 
The latter includes the optical selection rules (for details 
see Ref. @) • The vector notation is essential to account for 



the polarization degree-of-freedom which allows to create 
a non-equilibrium pseudospin distribution due to optical 
orientation [9(. 

The Hamiltonian describing the phonons and the 
carrier-phonon interaction is given by 
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FIG. 2: Sketch of constant energy contours of an electron 
subband (dashed lines), spin-split due to Rashba spin-orbit 
interaction. The dependence of the pseudospin orientation 
on the wave vector is visualized by the arrows. 



Kphonon = X^ W ( q ) &t ( q ) & ( q ) 

q 

+ E j E {ak m« (q) <L» ( k ' + q ) fc ( q ) c ^ c qo + h, 

+ E (.9^ m ;'(q)^(k' + q)6(q)^(k') + h.c.) \ , (4) 



with the annihilation (creation) operator of a phonon 
b") (q) . The linear interaction of phonons with elec- 
trons (holes) is ruled by the matrix elements g^, m »(q) 

{.9m 1 m"(q)) ■• wmcrl due to our choice of the energy eigen- 
states as basis is labeled by pseudospin indices. This 
expresses the fact, that a change of the wave vector 
due to a scattering event is in general accompanied by 
a change of the pseudospin state as visualized in Fig. 
|3J We note that without spin-orbit coupling the eigen- 
states are simple products of orbital and spin (-up or 
-down) eigenstates. Consequently, the matrix elements 
of the (spin-conserving) electron-phonon interaction re- 
duce to g m ' r m ' c ' = g S m > a m » . The actual dependence of 
the interaction matrix element on q is determined by the 
interaction mechanism, which has to be specified for a 
quantitative analysis [flLl- We note in pas sing, that the 
same system has been studied recently |21| even includ- 
ing besides electron-phonon interaction also impurity and 
carrier-carrier scattering, but not using the eigenstates 
of Ho . Instead the spin-orbit coupling is treated there 
explicitly as an effective k-dependent magnetic field giv- 
ing rise to an inhomogeneous broadening. In Sec. Ivl we 
shall discuss this situation, which is related to the present 
choice of basis by a k-dependent unitary transformation. 

In order to achieve expressions for the carrier-phonon 
scattering, we have to evaluate the EOM of the density 
matrix. In the basis of energy eigenfunctions of TLq the 
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FIG. 3: Possible scattering processes on the energy contours 
(thin full lines) connected with a change of the pseudospin. 

density matrix for the assumed 6-level model 

/ e ( mc m c ) (k) £(m cm „) (k) X 

falls into different blocks, where £»( m <=™ lc )(k) is the 2x2 
matrix for the lowest conduction band and ^( m « m »)(k) 
represents the 4x4 matrix for the hole states. The off- 
diagonal blocks g( m " m -)(k) and m <0 (k) include the 
coupling between valence and conduction band states 
by the exciting light field They describe the inter- 
band or optical coherence between the hole and elec- 
tron states coupled by the light field and their time 
derivative defines the EOM of the interband polariza- 



4 



tion P ma m „ (k) = (<4 lc (k) i4„ ( k )) • Without electron- 
phonon interaction the EOM of all entries of the 6x6 
density matrix would form a closed set of equations rep- 
resenting the coherent spin-dependent OBE for the sys- 
tem. A detailed theoretical study of the optical coher- 
ence and polarization dynamics, yet without addressing 
the spin/pseudospin, can be found in Ref. Il4l 

The pseudospin dynamics, in particular the relaxation 
and decoherence, is contained in the time evolution of 
the diagonal blocks, which shall be exemplified here for 
the electron system. The same steps of calculation would 
lead to the corresponding equations for the hole system, 



which however are more complicated due to the addi- 
tional orbital degrees-of-freedom. The 2x2 pseudospin- 
density matrix for the electrons is 

g(m c 7h c ) nA f Qm c rn a {^) Qm a — m c (k) \ 

V Q— m c m c (k) £>_ mc _ me (k) J 

The single entries are expectation values of products of 
a creation and an annihilation operator g mc fn c (k) = 
(4 (k) Cm c (k)) . We evaluate the commutators of the 
system Hamiltonian, Eq. (JIJ, with 4 C (k) Cm a (k) and 
take the thermal expectation value to find their EOM 



ihd t g„ lcf n c (k) = (e mc (k) - e mc (k)) p mcSlc (k) 

+ ]T {E(t) • d™ mv (k) p mc mv (k) - e* (*) • d-; m „ (k) pt c mv (k) } 

+ E m c (q)( c r< ( k + q) fo (q) ^ (k)) - <4 C , K (q) (4 C (k) &(q) c< (k - q)> 

q mj, 

+ <4* K (q) <4' c (k - q) 6 f (q) C^ c (k)) - I7lc (q) (4 e (k) (q) c m , (k + q)) j . 

(7) 

I 



The first two lines are the single-particle contributions 
of the SBE in Ref. @: they describe the dynamics 
caused by the spin-split energy levels and by the ex- 
citation of the electrons of either pseudospin from the 
valence subbands depending on the polarization of the 
driving light field. The three-operator terms specify the 
scattering of an electron (in one of the spin-split sub- 
bands) from one k to another one (in the same or the 



other spin-split subband) thereby absorbing or emitting 
a phonon, as visualized in Fig. [3] The three-operator 
terms (or their thermal expectation values) establish the 
phonon-assisted density matrix (|, whose entries obey 
EOMs of which we present as an example the one for 
s m > c m c (k + q, q) = {c] K (k + q) 6(q) c(k^J) 



ihd t s m > c rn c (k + q, q) = (e m / (k + q) - e^ c (k) - tiw(q}) s m / m a (k + q, q) 

+ E E {ah < ^ c k ( k + * + <*') b ^ b ^ c ^ W> 

k'q' 

- <4 C Ac (q')<4 4 (k + q) 6(q') &(q) ^ (k - q ')> 

+ g% ^ (q')(4 c (k + q - q') & f (q') %) ^(k)) 

- gfiZ ^ (q')(cl, (k + q) 6(q) &+ (q') c*, (k + q')) 

+ <4* fh c (q) (4. ( k <L' e (k + q) * t (k' + q) c„- lc (k)) 

-A> v rn v (q)(4. (k') 4, (k + q) V* (k' + q') c rTlc (k)) } . 

(8) 

I 



As can be seen from Eqs. Q and JSJ), we run into a increasing number of operators, which is typical for sys- 
hierarchy problem with EOMs containing terms with an 
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terns with interactions. This hierarchy problem can be 
overcome by a proper truncation. The standard proce- 
dure is to neglect the existence of coherent phonons cor- 
responding to the expectation value of a single bosonic 
operator (first order factorization) and to take into ac- 
count only the expectation values which lead to a phonon 
occupation number 0, Q • 

III. THE BOLTZMANN LIMIT 

The goal of the truncation is to gradually filter out the 
scattering terms up to a certain order in the interaction 
relevant for the investigated dynamics. To express the 
scattering contributions in the Boltzmann limit caused by 
electron-phonon interaction, we formulate the following 
rules for the truncation: 

1 . After factorization of the four-operator terms only 
expressions containing a macroscopic expectation 
value are taken into account. 

2. Scattering terms contributing in the Boltzmann 
limit are those proportional to the squared absolute 
value of the interaction matrix element in Eq. JTJl. 
This means that we neglect the so-called "polariza- 
tion scattering" due to inter- and intraband pro- 
cesses [2^ for which we refer to the next section. 

Applying these rules modifies Eq. (|SJ) and leads to an 
equation with the following characteristic structure 

d t x{t) =-iu)x{t)+y{t) , (9) 

I 



It has the characteristic form of expressions obtained 
from Fermi's Golden Rule: all terms are proportional 
to the absolute squared value of the interaction matrix 
elements and to the (^-function to warrant energy conser- 



where x(t) stands for the three-operator term and y(t) 
corresponds to products between phonon occupation 
functions /3(q) = (^(q) b(q)) and entries of the electron 
density matrix £> mc m c (k'). As presented in Ref. @, equa- 
tions of this type can formally be integrated to yield 



x(t) = :r(to)e~ Mt -* o) + / ° e~ luJt ' y(t - t') dt' . (10) 
Jo 



Inserting this result into the EOM of g nic m c (k) leads 
to a non-Markovian integro-differential equation, which 
can be solved analytically by applying the Markov and 
adiabatic limit [52, H3 ■ This corresponds to use instead 
of Eq. I|10|) the following expression 



x(t) = [-i-+«6(u))y(t), (11) 



where V denotes the principal value. With this approach 
we solve the EOM for the different entries of the phonon- 
assisted density matrices appearing in Eq. |7|. 

As we are interested in the relaxation due to electron- 
phonon interaction we write down only the scatter- 
ing contributions to the EOMs of the entries of the 
pseudospin-density matrix, Eq. ©, which correspond 
to damping terms. We obtain for the diagonal entries 



(12) 



(13) 

I 

vation in the scattering process. T™ cTOc (k) has the same 
form but with changed phonon and electron occupation 
factors. 



d t g mc mc (k)U = -IT* mc (k) g mc mc (k) + mc (k) (1 - g mc mc (k)) , 

I 

with the characteristic rates for scattering "out" of or The explicit form of r^ 4 OT (k) is given by 
"in" to the state with pseudospin m c and wavevector k. 

I 

Cfm c (k) = J Yl l^m c (q)| 2 { ( 5(e I< (k + q)-e mc (k)-^(q)) (l - £> m , m , (k + q)) /3(q) 

+ *(e mi (k-q)-e me (k)+Mq)) (l - ft„-<(k - q)) (l + /3(q))}. 
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For the off-diagonal entries we may write the scattering with 
contributions as 

d t Qm c -m e (k)| soo «l = _ r^ P _ me (k) Q mc -m c (k) , (14) 
I 



r^7- mc (k) = J E {b<m c (q)| 2 * 0<(k + q) - e-™ c (k) - ^(q)) X 

q m' c 

[(l - (?< m' c (k + q)) /3(q) + ffm' c m' c (k + q) (1 + /3(q))] 
+ l5m imc (q)| 2 <5 (e m ,(k- q) - e_ mc (k) + fkd{q)) x 
[(1 - Qm> c m' c (k - q)) (1 + /%)) + ^ m , (k - q) /3(q)] 

-m c (q)| 2 * (<W c (k + q) - £ ™ c (k) - fiw(q)) x 
[(l-e«(k + q)) /3(q) + ^m' c m' c (k + q) (1 + /%))] 
+ - mc (q)| 2 * (em' c (k - q) - e„ lc (k) + ^(q)) x 

[(l-0«(k-q)) (l + /?(q))+£> m , m ,(k-q)/?(q)] } . 

(15) 



Following the line of arguments in Ref . or 0, where 
the SBE have been derived for the two-level system of a 
conduction and a valence band state with carrier-carrier 
interaction, we may identify the damping rates in our 
pseudospin system with the inverse relaxation times 

J- = E ( r < < ( k ) + T ™U (k)) (16) 

J- = r e m 7_ mc (k) (17) 

of a Bloch vector, whose components are defined in the 
usual way by entries of the 2x2 pseudospin-density ma- 
trix. In spite of the similarities in the microscopic expres- 
sions of k and ?2,k for our system, we do not find the 
relation 2 Ti k = T2, k as for the system studied in Ref. 0. 
We notice instead a one-to-one correspondence between 
the individual contributions to both rates except for the 
sign changes in the pseudospin index m c , which - given 
the small spin-splitting - leads to a relation Tik — T^k • 
In fact for a system with spin-degenerate electron states, 
e„ lc (k) = e_ me (k) (as in systems with inversion symme- 
try), we find exactly Tik = T^k- This is in accordance 
with the argument used for two-level systems (see e.g. 
chapter 4 of Ref. E3), that a significant difference in 
these times (Ti 3> T2 ) arises, if due to separation of the 
two levels the decay of the population inversion requires 
energy dissipation. Making the gedanken experiment by 
assuming that the spin-orbit interaction is completely 
"switched-off" , the pseudospin and the orbital degree-of- 
freedom decouple. Hence, all dependencies on pseudospin 
indices are redundant as only spin-conserving scattering 
processes are possible. As a consequence no relaxation 



to a pseudospin-equilibrium is possible, because the scat- 
tering due to phonons is no longer capable of changing 
the pseudospin orientation. Nevertheless the scattering 
does not vanish but leads to a redistribution of the states 
in k-space within the separate pseudospin reservoirs. 



IV. BEYOND THE BOLTZMANN LIMIT 

In order to include all scattering terms up to second 
order in the electron-phonon interaction we have to go 
beyond the Boltzmann limit by taking into account also 
those contributions to the EOM of the phonon-assisted 
density matrices (taking Eq. JSJ as an example), which 
were omitted in the previous section. This is achieved by 
relaxing the second truncation rule used in Sec. IIIII and 
leads to additional contributions only to the EOM of the 
off-diagonal entry of the 2x2 electron pseudospin-density 
matrix which can be cast into the form 

<9 t (? mc - mc (k)| scatt 2 = -tt E £w o P -m' c (q) Bm' c -m' c (k+q) . 

(18) 

In contrast to Eq. I|14|) one has to sum here over the pseu- 
dospin index and the wave vector which enters differently 
in the self-energy £Q p „ m ^(q) and in Qm f c -m' e (k + q) • A 
corresponding scattering contribution was found in Ref. 
[H |(| for the interband polarization, i.e. for the off- 
diagonal entry of the 2x2 density matrix considered 
there. In order to present the structure of the self-energy 
we extract all contributions containing (according to Eq. 
HI 11 1) a S- function by writing 
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7T f 

f^_ K (q) = ^{5-m -mi(q)5mtmi(q) <J ( e -m' c (k + q)-e roc (k) + fia;(q)) x 

[(l- ftncmc (k)) /3(q) + g mcmc (k) (1 +/%))] 

+3-m c -m; (q) ffmt m£ (<j) 5 ( e < ( k + l) ~ £ -™c (k) + M<l)) X 
[(1 - p_ mc _ mc (k)) /3(q) + g_ mc _ mc (k) (1 + (3(q))] 

+9 e m' c mM) 9-*m' c -rajl) 5 ( e < ( k + <0 ~ £ -™<= ( k ) ~ M<0) X 
[(1 - g- mc -m c (k)) (1 + /%)) + e-™ c _ me (k) /3(q)] 

+ff<m c (q)9-<-m c (q) 5 ( e -mi(k + q) - £m c (k) - fiw(q)) X 

[(l-em eroe (k)) (l+/3(q)) + emcmc (k)/3(q)]} . 



We can identify the source terms composed of products of 
distribution functions and related to the different scatter- 
ing processes. As before, the energy conservation is con- 
tained in the delta function (therefore, we have omitted 
the contribution with g^, ^ (q) as a factor, because the 
energy difference between conduction and valence band 
states is usually much larger than the phonon energy). In 
contrast to Eq. I|15|) . the terms are not proportional to 
the absolute squared values of the interaction matrix ele- 
ments, which is typical for the contributions beyond the 

I 



dt Qra c m c (k) 



dt £>m c -m c (k) 



It contains in a microscopic formulation the pseu- 
dospin dynamics in electron subbands due to spin-orbit 
coupling, spin-selective optical excitation and electron- 
phonon interaction (for carrier-carrier interaction see the 
remark at the end of this paper) . By properly defining a 
Bloch vector as in Sec. IIHI and looking at the damping 
terms in the corresponding Bloch equations we can again 
specify the longitudinal and transverse pseudospin relax- 
ation times. As it turned out, only T2 : k is modified by 
additional terms (beyond the Boltzmann limit) discussed 
in this section, while 7\ k remains unchanged. 



(19) 

I 

Boltzmann limit as can be seen by comparing with the 
corresponding result for the two-level system of Ref. 6 
which shows a similar structure. There the terms beyond 
the Boltzmann limit have been denoted as polarization 
scattering [2^ with reference to the interband polariza- 
tion, while here they mean the corresponding scattering 
processes in the dynamics of g mc _ mc (k) . 

Together with the results of the previous sections we 
may now write the full set of EOMs for the pseudospin- 
density matrix 



(20) 



(21) 

I 

V. CHANGING THE SPIN BASIS 

When describing experiments designed to measure the 
spin relaxation time tsr and the spin decoherence time 
tsd of a system (see e.g. Ref. Il3land references therein), 
a basis is used with spin states oriented relative to a fixed 
direction, e.g. the growth direction of the QW structure. 
According to this choice, spins are spin-up (|) or spin- 
down (|) when aligned parallel or antiparallel to this di- 
rection, but in the presence of spin-orbit interaction spin 
is not a good quantum number. Consequently, the kinetic 
part of the Hamiltonian (including spin-orbit terms) for 
a general wave vector k is not diagonal. In order to be 



= Jn E i E W • d ^ m, (k) Pm c m v (k) - h.C.} 

-r£* ma (k) e mcmc (k) + r£_ mc (k) (l - g mcmc m 

= tt (e mo (k) - e_ mo (k)) £> mc _ mo (k) 
% n 

+ TE £ { E W ' rn v (k) Pm c m v (k) - E*(t) • d™* m „ (k) ^ (k)} 

(k) + ^ E mf-m' c (q) em' e -m' c (k + q) • 
qm' c 
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consistent with this convention, we translate the results 
of Sects. Illll and lTVI formulated in the eigenstates of Ho , 
to the spin-up/down basis. The unitary transformation 
connecting the two basis systems depends on the wave 
vector k and the type of spin-orbit interaction to be con- 
sidered. To keep the discussion as general as possible, 
we take into account the two most frequently discussed 
mechanisms of spin-orbit coupling, namely the linearized 
Dresselhaus term and the Rashba spin-orbit interaction 
pil |25| . Accordingly, we have instead of Ho the Hamil- 
tonian 

H] 1 — Hkin + Hr + Hd 1 (22) 

with the kinetic energy Hkin = k 2 ■ I2X2 of the 
free electron with effective mass m* . The Rashba- 
Hamiltonian has the form Hr = a (k x a y — k y a x ) with 
the Rashba coefficient a , the Pauli spin matrices o~ x / y 
and the components k x / y of the in-plane wave vector. 
The linearized Dresselhaus-Hamiltonian has a similar 
form, given by Hd = P (k x o~ x — k y o~ y ) with the weighting 



parameter (3 . The appearance of the Pauli spin matri- 
ces in Hr and Hd indicates the use of a basis with spin 
orientation parallel (or antiparallel) to the z-axis (which 
usually is the growth direction of the QW). The unitary 
transformation we are looking for is obtained by diago- 
nalizing H-\ 1 (Eq. J23) to find the eigenvectors 

with 



^{o?+(3 2 ){kl+kl)-±af3k x k y 

and the common abbreviation k± = k x + i k y . 

Applying the transformation matrix, composed of 
these eigenvectors, to the density matrix g( m ° m ^ (k) , we 
get the spin-density matrix in the basis of the spin- up and 
spin-down states 



,(U) 



(k) = 



£>n( k ) eu(k) 

£>!T( k ) <?ll( k ) 

d+(k) + 23?{A* g mc 
-d_(k) -2i3{A k g„ 



, c (k)} -d_(k)+2i3{^ mc _ mo (k)}\ 
. m „(k)} d+(k)-2 3?{A*g mc _ mc (k)} J 



(25) 



with d± (k) = g mc ma (k) ± Q- mc -m c s ) ■ For a particular 
choice of the spin-orbit interaction (Rashba or Dressel- 
haus) the corresponding unitary transformation can be 
derived on the basis of this result. 

Spin-dynamics experiments, such as time-resolved pho- 
toluminescence or Faraday rotation (for an overview of 
recent experiments usi ng t hese techniques see Ref. ITo|) 
or photogalvanic effect |26j| do not aim at the dynamics 
of the density matrix of an individual k but at quantities 
such as the spin polarization 

S = £( m (k)- m (k)) (26) 

k 

and the spin coherence |2l) 

C = 5> Ti (k)|, (27) 

k 

defined for the whole population of the two-level system. 
Their decay is characterized by the spin relaxation time 
tsr and the spin dccoherence time t$d ■ With the re- 
sults of Sects. IIIII and IIVI we are now in the state to 
formulate the relation between these quantities and the 
relaxation times ^ and T2 i k by applying the unitary 



transformation to express 

S = ^4K{A^ mo _ mo (k)} (28) 

k 

C = ^ v /d 2 _+43{A*p mc _ mc (k)}. (29) 

k 

The time-derivatives of S and C depend on those of the 
original pseudospin-density matrix g( mcmc )(k). Thus, 
the decay times of S and C , i.e. the spin relaxation 
and the spin decoherence are determined by T^k and 
T2,k derived in the previous sections, yet in a compli- 
cated relation. A microscopic calculation of tsr and tsd 
has to make use of this relation. Nevertheless, it is pos- 
sible to state, that the decay of the spin polarization 
S is determined only by Qm a —m, 

, (k) , i.e. by the trans- 
verse pseudospin relaxation time T2 j k, while the decay 
of the spin coherence C depends on both the longitudi- 
nal and transverse pseudospin relaxation times Ti,k and 
?2,k PH- The existence of a complicated relation be- 
tween Ti.k , 72.k for a simple spin-split two-level system 
and the spin polarization and spin coherence decay times 
of a whole carrier population has been mentioned before 
in the literature (see chapter 4 of Ref. [l(J but without 
making it explicit. 
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VI. CONCLUSIONS 

In this paper we have presented a microscopic formula- 
tion of spin dynamics in semiconductor heterostructures. 
It is based on the density matrix approach and its par- 
ticular form, the optical Bloch equations. Starting from 
the 6-level system of conduction and valence band states 
driven by optical excitation and including carrier-phonon 
interaction we derive explicitly the EOM for the 2x2 den- 
sity matrix of the electron subsystem whose energy lev- 
els are spin-split due to spin-orbit coupling. We employ 
a truncation scheme to include electron-phonon interac- 
tion in second order. In this limit we derive microscopic 
expressions for the longitudinal and transverse (pseudo-) 
spin relaxation times for the individual spin-split two- 
level system at a fixed k . Finally a connection between 
these results and spin relaxation times characterizing the 
dynamics of a whole population and accessible by exper- 
iments is established. It takes into account the differ- 
ent sets of eigenstates used in our microscopic derivation 



(which diagonalizes the spin-orbit coupling) and in the in- 
terpretation of the measurable times (with a fixed axis for 
spin quantization and nondiagonal spin-orbit coupling). 
Thus we provide at the same time a microscopic formu- 
lation of spin dynamics and its relation to experiments. 

We would like to emphasize that the concept presented 
here can be extended to include also carrier-carrier inter- 
action thus arriving at an extension of the coherent SBE 
of Ref. B For preliminary results we refer to Ref. I27I 
Further steps will be numerical evaluations of the micro- 
scopic expressions for realistic quantum structures and 
the explicit treatment of the spin dynamics for the hole 
system. 
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